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We establish bounds for the covariance of a large class of functions of infinite variance stable 
random variables, including unbounded functions such as the power function and the logarithm. 
These bounds involve measures of dependence between the stable variables, some of which are 
new. The bounds are also used to deduce the central limit theorem for unbounded functions of 
stable moving average time series. This result extends the earlier results of Tailen Hsing and the 
authors on central limit theorems for bounded functions of stable moving averages. It can be 
used to show asymptotic normality of wavelet-based estimators of the self-similarity parameter 
in fractional stable motions. 
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1. Introduction 

Consider a (non-Gaussian) symmetric a-stable {SaS, in short), a G (0,2), two- 
dimensional random vector (^,?7). By definition (see, e.g., Samorodnitsky and Taqqu 
[19]), this means that any linear combination uS^-\-vri is an SaS random variable, that is, 
it has a characteristic function £'c'^("?+'"') = q~'^"\s\" ^ where cr > is the scale coefficient, 
denoted by ||it^ -I- wTyHo,. We will sometimes assume that a > I so that the variables ^ and 
rj have finite mean. These random variables, however, always have infinite variance since 
a < 2. Nevertheless, functions K{(^) and L{rii) of these random variables may have finite 
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variance, as is the case for bounded functions K, L or unbounded functions K{x) = \x\^^^ 
with (3 < a/2, L{y) = In \y\. For such functions, the covariance 



is well defined. 

We are interested here in bounds on the covariance (1.1) when ^ and rj are dependent. 
(Examples of such ^,77 are provided in Section 2.) Wc seek bounds on (1.1) in terms of 
workable measures of dependence between ^ and rj. Since ^ and 77 have infinite variance, 
dependence measures other than the usual covariance must be used. Some alternative 
measures of dependence are available in the stable case, such as codifferencc and covari- 
ation (Samorodnitsky and Taqqu [19]). We will use a measure similar to the covariation 
and also introduce a new measure of dependence. 

Bounds on the covariance (1.1) are useful for several reasons. From a general perspec- 
tive, they complement the work on bounds of covariances of other random variables. 
In some cases, for example, Gaussian variables and integer power functions, these co- 
variances can be computed explicitly in terms of covariance between the two variables. 
In other cases, more general dependence measures, such as various types of mixing, are 
studied. As with bounds on covariances for other variables, bounds on (1.1) are useful for 
establishing limit results. In fact, we will use the obtained bounds to deduce the central 
limit theorem (CLT) for partial sums 



where K is as in (1.1) and is an SaS moving average time series with a € (1,2). As 
K can be unbounded, this result extends the CLTs for partial sums (1.2), with bounded 
functions K considered in Hsing [11] and Pipiras and Taqqu [16]. A related paper of 
Wu [23], which allows functions K to be unbounded, will be compared with our result 
below. We also apply the obtained CLT to establish the asymptotic normality of wavelet- 
based estimators of the self-similarity parameter in linear fractional stable motion. The 
obtained asymptotic normality result is, to our knowledge, the first of its type. 

The paper is organized as follows. Section 2 contains some preliminaries and statements 
of the main results on bounds of the covariance (1.1). A useful measure of dependence is 
studied in Section 3. The proofs of the main results on covariance bounds can be found in 
Section 4. Section 5 contains auxiliary results. In Section 6, we show the CLT for partial 
sums (1.2) and in Section 7, we apply it to deduce asymptotic normality of wavelet-based 
estimators in linear fractional stable motion. 



Cov(A'(C), L(r;)) 



(1.1) 




(1.2) 



n=l 



2. Results on covariance bounds 



Let (^,77) be an SaS, a g (0,2), random vector, as in Section 1. We will be concerned 
throughout only with distributional properties of {£.,r]). Therefore (Samorodnitsky and 
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Taqqu [19]), we may suppose, without loss of generality, that 

i^,v) I /(s)A/(ds), / 5(5)M(ds)Y (2.1) 



where (S*, /i) is a measure space, /, g S L°'{S,fi) and M is the so-called SaS ran- 
dom measure on S with the control measure ^. Heuristically, M(ds), s € 5, can be 
viewed as a sequence of independent SaS random variables with scale coefficients 
^(ds). More precisely, for any Borel set A, M{A) is an SaS random variable with 
scale coefficient ||M(yl)||Q. = /^(A)^/". Moreover, if Ai,...,An are disjoint Borcl sets, 
M{Ai U • • • U An) ^M{Ai)^ h Af (A„) and M(Ai), . . . , M{An) are independent. Be- 
cause the random variables M{Ai), . . . , M{An) are SaS, the resulting vector (^,77) is 
SaS. (For more details, see Chapter 3 in Samorodnitsky and Taqqu [19].) With the no- 
tation ||/||„ = (/g|/(s)|"Ai(ds))i/", the relation (2.1) also means that the characteristic 
function of (^, 77) is given by 

that is, -1- vriWa — \\uf + vg\\a- 

According to Samorodnitsky and Taqqu [19], the space {S,fi) in (2.1) can be taken as 
S* = (0, 1), /i = Lebcsgue measure. Other choices are possible and are often used. 

Example 2. 1 . Consider SaS random variables 

00 00 

k— — oc /c— — 00 

where ^ |/fc|" < 00, ^ |<?fc|" < 00 and {£„} is an i.i.d. sequence of SaS random variables 
(with scaling coefficient 1). Then, S — 1^, fi = counting measure, f{n) — fn and g{n) = gn- 

Example 2.2. Consider an SaS stationary time series given by the moving average 
representation 



= j ain-x)Midx), (2.3) 

where a G i"(R, da;) and M is an SaS random measure on R having the Lebesgue control 
dx. Take = ■^o and 77 = Then, (5, /i) = (R, dec) and /(•) = a{—-),g{-) = a(n — ■). 

We shall use the following quantities related to ^ and rj. 

Definition 2.1. Set 

Js 

[e,77]2= / |/(s)5(s)|"/V(ds). (2.5) 
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Observe that, unlike [Ci??]!! the quantities [C,77]i and [^,'7]2 are symmetric in their 
arguments ^ and ry. Observe also that, by Holder's incquahty, [^,?7]2 < oo with a G (0,2) 
and 77]^^ < cxD with a G (1, 2). The quantities r]]i and ri\2 can be viewed as measures 
of dependence between the random variables ^ and 77. In particular, the dependence 
measure [S,, r]]i is related to the covariation measure of a stable random vector (^, 77). The 
covariation, however, does not involve absolute values (see Samorodnitsky and Taqqu 
[19], Section 2.7). Observe, also, that the representation (2.1) of random variables ^ and 
r] is not unique: for example, another representation could be obtained from (2.1) by a 
change of variables. It is therefore important to note that the quantities in (2.4)-(2.5) are 
invariant with respect to the representation (2.1) of ^ and 77, as stated in the following 
proposition. See Section 4 for its proof. 

Proposition 2.1. The quantities [£,,11]*, '^'^'^ [^^vh in (2.4)-(2.5) do not depend 

on the representation (2.1) of{S,,r]). 

We shall often assume that ^ and 77 satisfy the following conditions. 

Assumption (Al ). There is a constant £1 > such that 

wfrJ^hrJ^- 1 i/(^)g(^)r/V(ds)>eiii/iis/2||gii2/2. (2.6) 

Assumption (A2). There is a constant £2 > such that, for all w,u G M, 

\\uf + vg\\l>e2{\\uf\\l + \\vg\0. (2.7) 

For example, the variables ^ and 77 in Example 2.2 satisfy these assumptions for suffi- 
ciently large n, where the constants £1 and £2 do not depend on n. Indeed, as proved in 
Section 4, wc have the following. 

Proposition 2.2. Let a G (0,2) and {^n}ngz be an SaS average time series from Ex- 
ample 2.2. Then, for sufficiently large n, the random variables C = Co o.nd r] = satisfy 
Assumptions (Al) and (A2), where the constants £1 and £2 do not depend on n. 

We shall now state the assumptions on the functions K,L in (1.1) which are used in 
this work. We shall assume the following. 

Assumption (Kl). There are f3 G (0,q;/2) and > such that, for \x\ > xq, 

|K(x)| <const|a;|''. (2.8) 

Assumption (K2). The function K{x)/x is non-increasing for x > Xq and non- 
decreasing for X < —xo ■ 
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Assumption (K3). The function K is bounded on \x\ < xq except possibly at a finite 
number of points xi^. . . ,Xp, around which the function K is integrable. 

Assumption (Kl) is quite natural since it ensures that EK{^i)^ < oo. Assumption (K2) 
is rather technical and depends on our method of proof. Examples of functions K satis- 
fying (K1)-(K3) are \xf with /3 £ (-1, a/2), (In |a;|)'" with m > 1 and many others. The 
choice of assumptions is also motivated by the applications to central limit theorems in 
Sections 6 and 7. 

The following results, proved in Section 4, provide bounds for the covariance 
Cov{K{^), L{r])) under various assumptions on the functions K and L. In Theorem 
2.1, we suppose that K and L are any two integrable functions. Theorem 2.2 concerns 
the functions K and L that are non-zero only for large values of the argument and which 
satisfy Assumption (K2). Theorem 2.3 is a combination of the previous two results. 
Finally, Theorem 2.4 is a consequence of Theorems 2.1, 2.2 and 2.3. 

Theorem 2.1. Let a G (0,2) and ^,77 be two SaS random variables (2.1) that satisfy 
Assumption (Al). Let also K,L G L"'^(R). Then, 

\Co^{K{0,L{r,))\<C\\KULUtri]2, (2.9) 

where the constant C depends on £1, a, and \\r]\\a- 

Theorem 2.2. Let a£ (1,2) and ^,7] be two SaS random variables (2.1) that satisfy 
Assumptions (Al) and (A2). Also, let K and L be two functions satisfying Assumptions 
(Kl) and (K2) with xq = b > \. Set Kb{x) = /ir(2;)l{|j.|>(,} and Lb{x) = L{x)\ ^\x\>b] : 
x^R. Then 

Cov(A'b(0,L6(77))| <C62^-"(K,??]i + K,^]2), (2.10) 
where the constant C depends on £1, £2, a, (3, ||^||q and ||?/||q. 

Theorem 2.3. Let a G (1,2) and ^,r] be two SaS random variables (2.1) that satisfy 
Assumptions (Al) and (A2). Also, let Kh be a function as in Theorem 2.2 and L G i^(M). 
Then, 

|Cov(Xb(e),L(r;))|<C6''-ilog(6+l)||L||i([77,e]^ + [77,e]2), (2.11) 
where the constant C depends on £1, £2, a, f3, \\^\\a and \\ri\\a. 

Theorem 2.4. Let a G (1,2) and ^,77 be two SaS random variables (2.1) that satisfy 
Assumptions (Al) and (A2). Also, let K and L be two functions that satisfy Assumptions 
(K1)-(K3). Then 

|Cov(X(0,i(r7))| < C([C,77]i + K,77]2), (2.12) 
where the constant C depends on £1, £2, K, L, \\£,\\a and \\ri\\a. 
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Remark 2.1. Observe that the bounds (2.10) and (2.11) hivolve the cut-off parameter 
b. In particular, as & — > oo, both bounds converge to 0, which is consistent with the fact 
that the two covariances also tend to 0. These results arc therefore acceptable in the sense 
that the effects of the cut-off b and the dependence are separated in the bounds. We also 
note that explicit expressions for the constants C in the bounds of the covariances above 
can be deduced from the proofs of the results, but these are not pursued here. 

3. Measure of dependence 

In this section, we establish some results on a measure of dependence U^^^ri between ^ and 
rj defined below. These results will be used in Section 4. 

Definition 3.1. For SaS random variables ^ and r/ in (2.1) and w,u gR, set 

U^^,,{u,v) = i;e"'«+''"' - Ee'^'^Ee'"' = c-ll"^+"f US _ g-ll^/HS-IK'sllS . (3.1) 

We can find the measure of dependence U^^ri in Section 4.7 of Samorodnitsky and 
Taqqu [19] , as well as 

lU^,v) = \\uf + vg\\2 - \\uf\\2 - \\vg\\2 (3.2) 

and — /^,,,(1, —1), which is called the codifference. Observe that 

\U^.jj{u,v)\<m,r,{u,v)\, u,veR, (3.3) 

due to the inequality |e~^ — e~^| < |a; — y\, x,y > 0. The next result, borrowed from 
Delbeke and Segers [10], provides various bounds on U^,r]- For notational simplicity, we 
shall write J/id/i for the integral h{s)fi{ds) below, whenever this is not confusing. 

Lemma 3.1. Let a G (0,2) and the function U^^rj be defined by (3.1). Then, for all 

\U^,r,{u,v)\<2\uvr'^[^,Tf^2, (3.4) 
\UU^,v)\ < 2|uz;r/V(l"l°'^ll«IIS^^-l'"l"'^ll''IIS^^)^[e,^]2, (3.5) 
\U^,^iu,v)\ < 2|uz;|"/V2(ll«IIS''ll''IIS''-[«'"]^)l™l°''[^,,7]2. (3.6) 

Proof. Inequality (3.4) follows from (3.3) and 

\lU^,v)\<2\uvr/' J |/5|"/MM = 2|^«|"/2[^,r7]2, (3.7) 

which follows from relation (5.8) below. As for inequalities (3.5) and (3.6), by using 
- 1| = |c^ -c"| <cl^^l|a;|, a; G R, and (3.7), we obtain 

|[/j_^(u,w)| =e-ll"-^IIS-ll''9llS|e-^«.''("''') -1| <e-""^IIS-ll"sll»|/^,,,(w,v)|e'-f«-'("'^)l 
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<2|uw|"/2e-(l"r^'ll/IIS^'-l''r''ll9llS^')'e-'(ll-^ll°^'ll»ll°^'-/ '^»l"^''^^)l™l°^'[^,77]2. 

This yields both inequahties (3.5) and (3.6) since /|/g|"/^d/x< ||/||q^^||5||q^^ and the 
exponents of the two exponentials are negative. □ 

The following two results concern the partial derivatives of the function U^^j]- We shall 
use the notation a^^^ = sign(a)|a|^ with a,p E M. The next lemma follows from Lemma 5.1 
below. 

Lemma 3.2. Let a G (1,2) and SaS random variables ^,77 be given by (2.1). Then, for 



(u,w) = -a (M/ + w.g)<"-i>/dAie-ll"-^+^fll = 



du 

+ a / (u/)<"-i>/dAie-ll"^IIS-ll''9llc 



(3.8) 



and 



^{u,v)^-a{a-l) / \uf + vgr-^fgAtic-\\^f+''s\\ 



du dv 

+ I {uf + vg)^°'-'^^fdfi y (w/ + i;g)<"-i>5d^c-""-'^+"f"° (3.9) 
'"(u/)<"-i>/dAi / (w5)<""^>5d^e-""^'llS-ll''9llo, 



a' 



provided, for the relation (3.9), that 

J\uf + vgr-^\f\\g\dfi<^. (3.10) 

Remark 3.1. When s G S* is such that /(s) = or g{s) = 0, we assume in (3.10) and 
(3.9) that \uf{s) + vg{s)\"~'^\f{s)\\g{s)\ = 0. This expression is clearly zero if one of /(s) 
or 17(5) is zero; if both are zero, then their values do not contribute in the definition (2.1) 
of ^ and ry. 

In the following lemma, we provide bounds for the partial derivatives in Lemma 3.2, 
using Assumption (A2). 

Lemma 3.3. Let a G (1, 2) and ^,77 be two SaS random variables given by (2.1). Under 
Assumption (A2), 



< C(|i;|"'i[r7,C]5:e-'^^(ll"«ll°+ll^''ll°) + \uf~^\U^^^{u,v)\), (3.11) 
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-iu,v) 



< C(C/i(u,«) + C/2(u,«) + U3iu,v)), (3.12) 



du dv 

where the constant C may depend on a,S2, \\C\\a o,nd W'qWa, o,nd 

Uiiu,v) = J \\uf\ - |«g|r-2|/||g|d/ie-'^(ll"«IIS+lt"'IIS), (3.13) 

U2{U, V) = (|up"-2 ^ |„|2a-2)g-e2(||<||S + ll«r/IIS) [^^ (3.I4) 

U3{u,v) = |ur->r-i|C/5,^(M,z;)|. (3.15) 

Proof. We consider only the bound (3.12), which is more difficult to prove. We shall 
denote by C a generic constant which may change from occurrence to occurrence and 
also by Ui, U2 and U3 the three terms on the right-hand side of (3.9). Since a — 2 < 0, 
by using \uf + vg\ > \ \uf\ — |wg|| and Assumption (A2), wc obtain 

\Ui\<C J ||^./|-|z;5ir-Vll3|dMe-^^(""^"^ + ll"«ll^)=CC/i(7.,^;). 

To bound \U2 + 11^1, add and subtract to U2 a similar term, where the first integral is 
replaced by J(u/)^"~^^/d/i, and to a similar term, where the exponential is replaced 
by e~""'^'*'"^l^'° . Then, by twice applying the triangle inequality, 

\U2 + U3\<C y|ur-i|/rdM J\v\''-^\g\''dfi\U^,„{u,v)\ 

+ C J \uf + vgr-'\g\dfie-\\^f+-<'\\' J \iuf + vgY'^-'^ ^ iuf)^"-'^\\f\d^l 

+ C J |7.|"-VrdA.G-ll"/+"fllS J \{uf + vg)<'^-'^ -{v9)''^-^^\g\d^i. 

Now, by using inequality (5.7) below and Holder's inequality, we obtain 

ic/2+c/3i<cii/ii^ii.9ii^i«ri|«r^ic/«,„(«,^;)i 

+ C(H"-i||/||r'|lg|U + |i;r-i||g||S)c-IW+''«IIS|„|-iy' \gr^\f\df^ 

+c\ur'\\me-\\"f+^^\\"«\ur' J \f^'\g\d^^ 

<C'{U3{U,V) + U2{U,V)), 

using Assumption (A2) and the relation < C(|up""2 + \v\'^"^^). □ 

In the following three results, we give bounds on various integrals of U^^jj and integrals 
of its partial derivatives. 
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Lemma 3.4. Let a G (0,2) and ^,77 be SaS random variables that satisfy Assump- 
tion (Al). Then 



(3.16) 



\U^^^{u,v)\dudv < C[^,r]]2, 
where constant C depends on Si,a, \\S,\\a o.nd W'qWa- 



Proof. We consider the integral over (0, 00) x (0, 00) only and examine it over four 
regions. Over (0, 1) x (0, 1), by using (3.4), we have 



1 pi 



Jo 



1 fi 



\U^^riiu,v)\dudv<2 / u°'/^v°'/'^dudv[tvh^C[tTi]2- 
Jo Jo 



Over (1, 00) X (1, 00), by using (3.6) and Assumption (Al), we can boimd the integral by 



oc />oo 



u^/^v^^'^dudvl^^rij^ 



up 



Observe that the integral here is finite, since we can bound e" 
to a constant by {uv)~p for arbitrarily large p> Q. Over (0, 1) x (1, 00), by using (3.5), 
we have as a bound. 



("''''\\f\rj'--''hrj'ru"/'v"/^dudv[^,7^h 



The integral here is again finite. A similar bound holds over the region (l,oo) x (0,1) 
and hence the result (3.16) is valid. □ 

Lemma 3.5. Let a€ (1,2) and ^,ri be SaS random variables that satisfy Assumptions 
(Al) and (A2). Also, let 



F{u) 



bf>-^\u\-\ 



if \u\ < 1/6, 
if \u\ > 1/6, 



(3.17) 



where b>l and (3 G (0,q;/2). Then, 



du 



{u,v) 



F{u) dudv < Cb"-' log(6 + mv,^]l + vh), (3.18) 



du dv 



{u,v) 



F{u)F{v)dudv < C62'3-"([^,r;]i + [^,7712), 



(3.19) 



where the constant C depends on a, (3, 61,62, ||CIU '^'^'^ !I'?IIq- -^^ particular, the derivative 
dvdu ^i'V defined a.e. ditd^;. 
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Proof. Consider first inequality (3.19), which is more difficult to prove. By Lemma 3.3, 
it is enough to show that the bound (3.19) holds when gf q,„ U^,n is replaced by Ui, by 
U2 and by C/3. Since Ui, U2, U3 and F are even functions in u and v, it is enough to 
consider the integral over (0, 00) x (0, 00). Let us denote this integral with Ui, U2 and C/3 
by /i, I2 and I3, respectively. That is, Ik = Uk{u,v)F{u)F{v)dudv, k = 1,2,3. 

As in the proof of Lemma 3.3, we shall denote a generic constant by C . 
Bounding Ii. We have 



00 /"OO 



h< / \f{s)\\g{s)\ 



\\f{s)\u-\g{s)\v\°' ^F{u)F{v)Audv^i{ds) 



= :/ \f{s)\\g{s)\h{s)^i{ds). 

JS 



By using Lemma 5.2 below, when /(s) 7^ and g{s) ^ 0, we obtain 

a-2 

F[u)F{v)dudv\g{s)\°'-'^ 



noo 

< Cb'^-^' { 1 



5(s) 



I{s) 



a-2 



9{s) J 

= c62''-"(i/(5)r-2 + |g(s)r2). 

In view of the remark following Lemma 3.2, we conclude that h < Cb^'^-'^[^,r]]i. 

Bounding I2: We consider I2 over four regions (0,1/6) x (0,1/6), (1/6, 00) x (1/6, 00), 
(0, 1/6) X (1/6, 00) and (1/fe, 00) x (0, 1/fe), and denote the corresponding integrals by /2,i, 
-^2,2, ^2,3 and I2A- Wc shall also consider U2 with only the term |w|2"~2 (^^^ ^^ie case of 
|?;p"~2^ the bound is obtained by symmetry). Over (0,1/6) x (0,1/6), after the changes 
of variables w w/6, u — > v/b, we obtain 

.1/6 i^l/b 



< 



since -/? + 1 > and 2a - 2 - /? + 1 > 0. Over (1/6, 00) x (1/6, c»). 



00 poo 



h.2 = 



y2Q-2g-£2i«ni/iis-e2ifnisii; 



3-l,,-U/3-l..-l, 



n/b Jl/b 
< 62'^-2 



^2a-3^-e,\ur\\nZ^^ 







1/6 



^b'^-'^v-^dudv[C,ri]i 
-e2knisllS„-id„[^,,7]i 



<cb^^-^iog{b+m,v]i 

since 2a - 2 > 0. Over (0, 1/6) x (1/6, 00), 

f.l/6 /.oo 



l2,3 < 



^^2a-2^-e,\vn\gK^- 



1/b 



^b^-\-^dudv[^,7]]i 
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^0 Jl/b 

<C6-2"+2/3iog(6+l)[e,r;]i. 
Similarly, over (1/6, oo) x (0,1/6), 

Jl/b Jo 

poo /"l/^ 

Jo Jo 

We conclude that h < Cb^>^-'^\og{b + l)[^,r]]i. 

Bounding I^: We shall here use inequalities (3.4)-(3.6). As in the case of I2, denote 
the integral I3 over the same four regions by I^^i, 13^2, I3.3 and 73.4, respectively. Then, 
by using (3.4), we have 

.i/b .1/6 

h.i<C / M""it;"-^W2„"/2u-'3v-'3dudw[e,?/]2<C6-3"+2f[^,,y]2. 
Jo Jo 

By using (3.6) and Assumption (Al), and making the changes of variables u = 
j^-i^2/a _ fj-iy2/a i,q\q^ ^ obtain 

/b 

noo 

By a further changes of variables x = b°'w and wy ~ z, we have 

/a 2 < C62/3-3a+a+2(a-l) r dw /'"e-2ei||/||S/^|ls|l°/^.^2(a-l)/aj^[^^^]^ 

Jb-" W Jw 
<CT2/3-2log(6+l)[^,^]2, 

by splitting and bounding < /o°° + Ii° ■ Turning to 13 3, this time us- 

ing (3.5) and making the changes of variables u ~ x^/°'b~'^,v ~ y'^/°'b~^ below, wc obtain 

rl/b poo 

h,3<CbP-' / / «"-iz;"-i«"/2^;"/2e-(«°''ll/ilS''-°''ll9llS^')\-^-idudi;K,r7]2 

^0 Jl/h 

Jo J I 

By making another change of variables yb~"/^ = w, it is easy to see that the last expres- 
sion can be bounded (up to a constant) by &2/3-3a^(a/2)-Ha-i ^ j,2/3-((3a)/2)-i[^^^]2_ 
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One obtains the same bound for 73.4 by symmetry. We can now conclude that < 

cb^^~Hog{b+m,v]2. 

The result (3.19) of the lemma follows from the bounds obtained for /i, I2 and I^. The 
inequality (3.18) can be shown in a similar way by using (3.11). □ 

Lemma 3.6. Let a £ (1,2) and ^,77 be SaS random variables (2.1) satisfying Assump- 
tions (Al) and (A2). Also, let 



G{u) 



\u\ 
1, 



where /3 S (0,q;/2). Then, 



du 



{u,v) 



if \u\ < 1, 
if \u\ > 1, 



G{u)dudv < C < 00, 



// 



{u,v) 



G{u)G{v) dudv < G < 00, 



du dv 

where the constant G depends on a, (3,Si,e2,\\^\\a arid ||?7||q. 



(3.20) 

(3.21) 
(3.22) 



Proof. We prove only the bound (3.22). The generic constants G and C" below are as 
in the statement of the lemma. As in the proof of Lemma 3.5, we consider the integral 
(3.22) over (0, 00) x (0, 00) only and denote it by Ji, J2 and J3 when gf g„ U^^r; is replaced 
by Ui, U2 and U3 in Lemma 3.3, respectively. To bound Ji, we can apply the bound for 
the integral Ii in the proof of Lemma 3.5 with 6=1 because the difference between the 
functions F and G can be accounted for by bounding the term e~^2'^ll"-''IIS + ll'"9lt°) in JJi by 
C|u|-i|t;|-i when |u| > 1 or \v\ > 1. We therefore obtain Ji < C[C, < C(||C||:^-i ||7?|U + 
lk/lla~"'^IICIIct) < C" < 00. In the case of J2, we have 

J2= / / (|u|2"-2 + |i;|2"-2^e~"^(""^"° + ll''f"S)G(M)G(w)dMdwK,r;]i 
Jq Jo 

and by bounding the exponential by C'e"^^'!!"'''!!"^!!'"^!!^) when |u|>lor|f|>l, 

we can again account for the difference between the functions F and G. We can therefore 
again apply the bound for the integral I2 in the proof of Lemma 3.3 with 6+1 and obtain 
J2 < C[^,'7]i < C" < 00, as above. The integral 

J3= / / \u\°'-^\v\°'-^\U^,n{u,v)\G{u)G{v)dudv 
Jo Jo 

can be bounded by using techniques from the proof of Lemma 3.4, together with the 
bounds (3.4)-(3.6) for and the inequality [^,r]]2 < ||^ l|a^^ ||r;||a^^ = C < cx). □ 



Bounds for the covariance of functions of infinite variance stable random variables 1103 

4. Proofs of the main results 



We shall prove here Propositions 2.1 and 2.2 and Theorems 2.1-2.4. 

Proof of Proposition 2.1. The vector (^,77) in (2.1) has the so-called minimal repre- 
sentation 

i,,,i{lMm^.ljim^^). (4.1, 

where M is an SaS random measure on S with the control measure Jl. For more infor- 
mation on minimal representations, see Rosihski [18] or [17], Section 2. It is therefore 
enough to show that the quantities [C,?7]|, [^,ri]i and [^,?7]2 in (2.4)-(2.5) remain the 
same when either the representation (2.1) is used or when the minimal representation 
(4.1) is used. We consider only the quantity [^,7?]2 in (2.5). 

We may assume, without loss of generality, that {s : f{s) ^ or g{s) ^ 0} = 5 ^-a.e. 
(otherwise, the set {s : f{s) = 0, g{s) = 0} = S can be eliminated from the representation 
(2.1) without changing t]]2)- By Remark 2.5 in Rosihski [17], there are maps (/>: S ^ S 
and /i : 5 1-^ R \ {0} such that 

fis)^his)ms)), ff(s) = Ms)5(0(s)) (4.2) 

a.e. /i(ds) and 

fl = ^ih o 4>~^ , (4.3) 

where /x/i(ds) = |/i(s)|"/i(ds). By using (4.2) and (4.3), and making a change of variables, 
we have 

\f{s)g{s)^/^^,{ds)= [ \Ms))g{cl>{s)r/'fi,{ds)= / |/(i)g(S)|"/2/I(d5). 
s Js Js 

This proves the result for the quantity [^,?y]2 in (2.5). □ 

Proof of Proposition 2.2. Let (T„a)(a;) = a{n — x), n £ Z, a; £ M. In view of (2.3), ^0 
and ^,1 have kernels / = Tqo and g = T„a. That ^0 and ^„ satisfy Assumption (Al) follows 
from the facts that ||r„a||Q = ||2oa||Q and |(Toa)(x)(r„a)(x)|"/^ dx as n — ^ cxo . To 
see why the last integral converges to 0, first write it as 

/•oo pi 

\a{-x)a{n - x)\°'/^ dx ^ / \a{~x)a{n - x)]""/^ dx + / \a{-x)a{n - x)\°'/^ dx 

J I J —00 



< 



ll«(-)l(.oo)li:/'l|a(-)||S/^ + ||a(-)IIS/'l|a(-)l(„-.oo)li:/', 



by the Cauchy-Schwarz inequality and a change of variables. Now, observe that the first 
term in the bound is arbitrarily small for large fixed I and that the second term converges 
to for fixed Z as n — > 00. 
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We now show that they also satisfy Assumption (A2). We consider only the case 
a e [1,2); the case a £ (0, 1) can be proven by considering || • jj^ below. For some I £ M, 
we have, by 2{x + yY >x^ + y^, x,y,p> 0, and by Minkowski's inequality applied to the 
norm || • \\a, 

\\uToa + vTna\\a = \\{uToa + vT„a)l(-ocj) + (wTba + wr„a)l(;^oo) L 

> l\\{uToa + vTna)l(^^oo,l)\\a + 5 11(^^700 + wT„a)l(;_oo)L 

— oc\l) Wa) 

^ (4.4) 

+ 2(IK'rnal(i,oo)L - \\uToal(i,^f\\J 

= sdl^oalc-oo.olL - \\Toal(i,oo)\\J\u\ 

— n,oo) 

|l - ||roal( 

)lla)l"l- 

Since ||Toa||Q — limm^oo ||Toal(-oo.m) lU? there is a constant M > such that, for all 
m>M, 

1 1 Tool ( )IL-||T^oal( 

m,oo 

)lla> 1^0^11 

1 1 Tool ( 

— m,oc 

By applying these inequalities to (4.4), by fixing I > M and taking n such that n — l>M 
or n > [M + /), we obtain 

\\una + vT^a\U > J(||^iToa|U + lK'T„a|U) > ^-l^(||^iToa||^ + |K,T„a||2)i/", 
which is Assumption (A2). □ 
We shall denote the Fourier transform and its inverse by 



(4.5) 



iTh)iu) = / e™-"/i(x)dx, 
{T"'k){x) = {gk){x) = -i- / e-'^-k{u) du 

and we shall also use the measure of dependence U^^rp which was analyzed in Section 3. 
The following two lemmas are used in the proofs of Theorems 2.1-2.4. 

Lemma 4.1. Let K he a function satisfying Assumptions (Kl) and (K2) with xa = h>l. 
Set Ki,(x) = K{x)\s^\j.\y^ anrf /<(,.„(x) = A'(x)l{b<|x|<n}; a;GR, nGN. Then, 
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where the constant C does not depend on b. Moreover, for fixed b, 



(u) 



<C 



if \u\ < l/b, 
ifueR, 



(4.7) 



where the constant C does not depend on n. 
Proof. The second inequality in (4.6) follows from 



27r 



,Kb{x) 



< 



< 



-fa 



K{x) 



4i\Kib)\ + \K{-b)\) ^ 
b\u\ 



where we used Assumptions (Kl) and (K2), together with the inequality 



e'^"/(x)da; 



< — max I /(x) I, 

\u\ xe[c,d] 



(4.8) 



valid for c < d, it G R and a fimction / monotonia on [c,d]. (The inequality (4.8) is a 
consequence of the second mean value theorem; see, e.g., in Kawata [12], page 24). As 
for the first inequality in (4.6), when < |u| < l/b, we obtain 



27r 



,Kbix) 



(u) 



<r I 1 1^(^)1 . 

< I l{&<|x|<i/|u|}— j— j— dx- 



K{z/\u\) 



g-izsign(M) J. 



Using the fact that "^^/l^"^"* is monotone for > b, hence for \z\ > 1, and apply- 

ing (4.8), we obtain 



27r 



(u) 



< 2 



l/\u\ 



K 



<C\u\ 



The bound (4.7) can be proven in a similar way. 



□ 



Lemma 4.2. Let a e (1,2) and ^,77 be two SaS random variables (2.1) that satisfy 
Assumptions (Al) and (A2). Then, for {x,y) gM^, x,y^Q, 



iGU,.)i^^y) = --[gg^^U,^,]ix,y). 



(4.9) 
(4.10) 



Proof. Since, by Lemma 3.6, -^U^^ri, g.^g^ U^^jj £ L {M. ) and also, by Assumption (A2) 
and the bound (3.11), lim„^oo U^^ri{u,v) = and lim„^oo(^J75^^)(it, w) = 0, we obtain, 
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using the integration by parts formula, that 



IX Jr du IX \ du J 



dudv □ 



We next provide the proofs of Theorems 2.1-2.4. 

Proof of Theorem 2.1. Let ip^^n, V'f and -0^ be the density functions of (C,'7), 6 a-nd 
77, respectively. Also, let 4>^.rj, 0{ and (pri be the characteristic functions of (^,?7), ^ and 
77, respectively. By using the inversion formula for the densities of (^,77), ^ and ry, and 
K,Le L^ilS.) and f7^^^ G L1(R2)^ we obtain that 

Cov(X(e),i(?7)) = EK{OL{rj) - EK{OEL{r^) 

^ K{x)L{y){ip^^r,{u,v) - -!/'«(w)^/'t7("))da;dy 

/ ir(x)L(y)((g</)5,,)(x,y)-(e05)(x)(e0,)(2/))d.Tdy 

i^(a:)L(y)(a[/c,„)(x,2/)da:dy 

(eA-)(w)(gL)(i;)C/^,^(ii,«)dud7;, 



R2 



by Fubini's theorem. The bound (2.9) follows by using \{gK){u)\ < C||ii:||i, \{gL){u)\ < 
C\\L\\i and Lemma 3.4. □ 

Proof of Theorem 2.2. For n > 1, set Kb,n{x) = Kb{x)l[\x\<n} = K{x)l{b<\x\<n} 
and Lb^nix) = L{x)l^b<\x\<n}-, ^ G II^- Since, by Lemma 3.6, q^q^ U^^^^ e i^(M^) and 
Kb^mLb.n g£^(K), we obtain, using (4.10), that 



Cov(ii:f,,„(^),L(,_„(77)) = / Kb^n{x)Lb,n{y){GU^^y^){x,y)Axdy 

I ^S^b^(g^u,^ix.y)dxdy (4.11) 
Js2 X y \ oudv J 
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y — ■ (m) y — ■ (^) a o iu,v)dudv. 

X J \ y J ouov 



Now, observe that, by using Assumption (K2) and the inequahty (4.8) above, we have 

1 ^^--^^(-)l,,<,,,dx=:(g^)(.) 



27r 

as n — > oo, for all u e M {u^ 0). Then, by letting n — > cx) in (4.11) and using (4.7) in 
Lemma 4.1, Lemma 3.6 and the dominated convergence theorem, we obtain that 

Kb{x)\, JrLb{y)\, ,d^U^„ 



CoHKbiOMm = - 1 (G^ ) in) ) d-d.. (4.13) 

The bound (2.10) follows from (4.6) in Lemma 4.1 and (3.19) in Lemma 3.5. □ 

Remark 4-1- Note that the proof of Theorem 2.1 fails for Theorem 2.2 because it may 
happen that QKf, is not defined (this is the case, e.g., when Kii{x) — log2 |a::|l||a;|>t,}). This 
is why we consider QKi,{x)/x instead in (4.12). Note, also, that {GKi,{x)/x)(u) in (4.12) 
is an improper integral because we could have Kbix)/x ^ L^(R) U L'^{R) and thus the 
usual L^- or L^-Fourier transforms for Kb{x)/x may not be defined. What helped was 
the monotonicity of the function Kb{x)/x. 

Remark 4. 2. In order to avoid dealing with complications involving g°g^, ^g.r;; one may 
think of transferring (by the integration by parts formula) the derivative q^q^ from C/j 
onto {Q-) to obtain 



CoY{Kb{O.Lb{ri))^- 



The problem with (4.14), however, is that the derivatives of {Q-) are not easy to manipu- 
late. Consider, for example, the functions K}j[x) = Lh{x) = |x|^l{|a;|>b} with (3 G (0,a/2) 
which appear in the wavelet applications. We can verify that 

d ( r^b{x)\ ■ I \( fD -0-1 /"°° 0-1 ■ J , ./jsinfeit 

— — C/ (m) = — sign(u) LOu ^ / smzdz + o*^ 

du\ X J n \ Jbu u 

and hence the right-hand side of (4.14) involves the integral with sin 6m sin 6u as part 
of its integrand. It turns out that we cannot estimate this integral by putting absolute 
values on the integrand (the sign of sinhusaibv is important). We therefore worked with 
the formula (4.13) instead. 
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Remark 4- 3- Although the formula (4.14) is not used here, it may be useful in other 
contexts. For example, if we take ln| • | for the functions Kf, and Lh, we obtain 

Cov(ln|e|,ln|7y|) = l / / U^^^v)^. (4.15) 
4 7rJr \u\\v\ 

The formula (4.15) is proved rigorously in Delbeke and Segers [10], where it is used to 
bound the covariance (4.15). 

Proof of Theorem 2.3. For ??. > 1, set Kb,nix) — K{x)l^i,<\x\<n}- x £M. Since, by 
Lemma 3.6, ^C/^.r; £ L^(R^) and Kh n G i^(R), we obtain, using Lemma 4.2, that 



Cov(if6.„(e),L(7,))- / Kt.n{x)L{y){gU^,^)ix,y)dxdy 

1 f Kb,n{x) J , n ^ 



1 ./b2 X 



L{y)(^g^U^,^{x,y)dxdy (4.16) 



1 / fn^bAx)\^^^^g^^^^^dU^ ^^ ^^^^^ 

ou 



By letting n ^ oo in (4.16) and using (4.12), (4.7), Lemma 3.6 and the dominated con- 
vergence theorem, we obtain that 

1 f fnKb{x)\,^^,,^^,,^^,dU^,„ 



Cov(/vb(0,L(??)) = T / [g-^]iu)igL)iv)^^iu,v)dudv. (4.17) 

The bound (2.11) follows from (4.6) in Lemma 4.1 and (3.18) in Lemma 3.5. □ 

Proof of Theorem 2.4. Let xo,k and xq^l be the points in Assumptions (K1)-(K3) 
corresponding to the functions K and L, respectively. For b > max{xo./<-, xq^l} and 6 > 1, 
set Ki,{x) = X(x)l{|j.|>ft| and Lb{x) = L{x)l^\x\>b}, a; G M. By Assumption (K3), we have 
X = iffc + X and i = Lf, + Z, with K,Le L\R, dx). The bound (2.12) follows by first 
writing 

\Cov{K{0,m)\ <\Cov{Kb{0,Lt{ii))\ + \Cov{Kb{OMv))\ 

^ ^ (4.18) 
+ |Cov(A-(e),ib(r,))| + I Cov(A'(e),L(r,))| 

and then using Theorems 2.1, 2.2 and 2.3 to bound the terms on the right-hand side 
of (4.18). □ 



5. Auxiliary results 

We prove here some results which were used in Sections 2 and 3. The following lemma 
was used in the proof of Lemma 3.2. It involves signed powers a^^^ = lal^* sign(a). We 
write J hd^ = /i(s)/i(ds) for notational simplicity. 
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Lemma 5.1. Let f , g £ L"" {S , iJ.) withal (1,2) and set 

F{u,v)= J \uf + vg\°'dfi, u,veR. (5.1) 

Then, for u, v S K, 



dF , 
du 



dv du 

provided, for the second equality in (5.2), that 



(5.2) 



(it,w) = a(a-l) / + '^fgA^l, 



I \uf + vgr-^\f\\g\dfi<c^. 



(5.3) 



Proof. By the mean value theorem, we have, for some u* = u*{u, v, h, x) € [u,u + h] (we 
consider h> 0), 



UF{^ + h,v)-F{u,v))~a f {uf + vg)^''-^^ f dpi 



^{\{u + h)f + vg\''-\uf + vgndfi^a j {uf + vg)^''-'^ f dfi 



{{u*f + z;g)<"-i> - {uf + «5)^""'^)/d/i 



<2ay \u-u*\°'-^\f\°'dpi<2ah°'-'^ J l/l^d^i^O, 



by (5.7) in Lemma 5.3, as /i — > 0, which proves the result for the first derivative. The case 
of the second partial derivative of F is more delicate. By the mean value theorem with 

V* = V* {u, v,h,x) £ [v,v + h], we have 

l(^^{u,v + h)-^{u,v)^^aj^{{uf + {v + h)g) - {uf + vg) <""!>)/ d/. 

= a[a-l) j \uf + v*gr~^fgdtx. 

Since \uf + v*g\°'~'^ — ^ \uf + w.g|"~^, as ft- J, 0, and. by (5.6) in Lemma 5.3 below, 
{a-l)\uf + v*gr-^\f\\g\ 
-\\{nf+{v + h)gf'-'^ - {uf + vg)^^-'^\\f\ < 2\uf + vgr'\f\\g\, 
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the dominated convergence theorem and condition (5.3) together imply that 

^(^^{u,v + h)-^{u,v)^!^a{a-l) J \uf + vgr' fgd^, 

which proves the resuh for the second partial derivative. □ 

The next lemma was used in the proof of Lemma 3.5. 

Lemma 5.2. Let a e (1,2), r > and F{u) ^ {ul-l^ for \u\ < 1/b and F{u) = b'^-^\u\-'^ 
for \u\ > 1/b, where b > I and f3 £ {0, a/2). Then, 

Io{r):^ / \ru-v\°'-^F{u)F{v)dudv<Cb^'^-"{l + r°'-^), (5.4) 
Jo Jo 

where the constant C depends only on a and (3. 

Proof. Observe that F{b~^u) — b^Fi{u), where Fi{u) — \u\~^ for |m| < 1 and Fi{u) = 
\u\~^ for \u\ > 1. Then, since 



/>oo />OC 

/ / \rb-^u~b-'^v\°'-'^F{b-'^u)F{b-^v)b-'^dudv 
Jo Jo 

/>oo />oc 

fe^'^-" / / \ru-v\"-'^Fi{u)Fi{v)dudv, 



it is enough to show (5.4) with b—1. Consider the integral Io{r) over four regions (0, 1) x 
(0,1), (1,cxd) X (l,oo), (0,1) X (l,cxo) and (l,cx3) x (0,1), and denote the corresponding 
integrals by /o,i(r), /o,2(''), ^o,3('') and /o,4(r), respectively. By making the changes of 
variables v = ry, u = yw below, we obtain that 

Io.i{r)= I I \ru-v\'^-^u-l^v-l^ dudv 



JO 
l/r / 1.1/y 







dy\ / dw\w -l\°'-^w-^ ]y 



First, suppose that r > 1. We examine /o.i('') hi three cases, depending on the behavior 
of the inner integral above. 

Case 1: a — 2 — /3 + 1 = Q! — 1 — /3>0. In this case, since r > 1 and a < 2, we have 

l-l/r 

h,l{r)<C / 2/-("-l-'3)y"-l-2/3^^^a-l-/3 
"'0 

Case 2: a — 1 — /3<0. In this case, since r > 1 and /3 < 1, we have 
loAr) < C f'^\''-'-^^dyr"-'-^ ^C'r^-' < C . 
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Case 3: a — 1 — /3 = 0. By using the integration by parts formula, we have 
hAr)<C / \\ogy\y-'^ dy<C'{r'^-^\ogr + rf^-^) <C" . 



This shows that /o.i('') < C when r > 1. If < r < 1, then, by writing /o.i(r) = 
lo So I" ~ r~^v\°'^^u~^v~^ dudvr°'~^ and using symmetry, we can conclude, as above, 
that /o,i(r) < Cr"-2_ Hence, /o,i(r) < C(l + r"'^^) for any r > 0. 

Turning to the integral Iq^t), supposing first that r > 1 and making the change of 
variables ru = y below, we obtain that 



\ru-v\°'-^u-\r^dudv = dy J dv\y - v\°'-^y-'^v-'^ < C. 



By using the symmetry argument above (see the case of /o,i(r)), we can deduce that 
loM^) < r"''^ when < r < 1. Hence, /o,2(?') < C(l +r"-2) for any r > 0. 
By making the changes of variables v = ry, u = yw, we obtain that 



= I I dv\ru — vl" u ' v 

,1/y 



'1/r 

Consider first r > 1. In this case, we have 



dy^ 'du;|u;-l|"-2u;-/5^y"-2-/3^a-2 



(5.5) 



IoAr)=J dyl^J du;|u;-ir-2w;-'5jy"-2-/5 r"-2 

^2^(/'^'du;|u;-ir-2u;-'')j/"-2-^ 3M + 3M- 

For the integral I^^ir), we have ll^ir) = C y/3-iy"-2-/3 dy r"-^ < C'r"-^ < C" since 
r > 1. To bound the integral /q 3(r), we consider the following three cases. 
Case 1: a — 1— /3>0. In this case, we have 

/o3W<C C y-("~i-^)y"-^-'^d2/r"-2^C"(logr)r"-2 <c". 

Jl/r 

Case 2: a — 1— /3<0. In this case, we have 

llz{r)<C r y''-^-Pdyr'--^=C'r^-^<C'. 

Jl/r 

Case 3: a — 1— /3 = 0. In this case, we have 

ll^{r)<C f |logy|y-Myr"-2 = c'(logr)V"-2<C". 

Jl/r 
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This shows that I^^ir) < C when r >1. Now, consider < r < 1. Since r^^ > 1, by 
using (5.5), we obtain that 

Hence, /o,3(r) < C(l + r"-^) for any r > 0. It can be similarly shown that the same 
boimd holds for /o,4(r). This proves the result (5.4). □ 

The following inequalities were used in the proofs of Lemmas 3.1, 3.3 and 5.1. 

Lemma 5.3. Let xi,X2 G R. When a G (1,2), the following inequalities hold 

\x["-'^ -Xi"~'^\<2\x2\"-^\XI~X2\ {X2^0), (5.6) 

\x["-^'> -xi''-''>\<2\xi-X2r-\ (5.7) 
Moreover, when a G (0,2), 

\\x + yr~\x\"-\y\''\<2\xyr/\ (5.8) 

Proof. The inequality (5.6) follows from the inequality |(1 + z)^''^ ^ 1| < 2|z|, with q S 
(0,1) and 2 e M, in Lemma 3.4 of Cioczek-Georges and Taqqu [7]. We provide here a 
geometric proof, which is best understood by graphing the function p{x) = x£R. 
Suppose, without loss of generality, that X2 < 0. Then, there exists an x* such that the 
slope of the line connecting the points {x2,x^^^'^) and (a;*, (.t*)^"^^^) is the largest 
among the slopes connecting {x2,x^ ^'') and {x,x^"~^^) with a: G M, x ^ X2. This x* 
satisfies < x* < |a;2|. Then, 

i4°-^>-4°-^>i ^ i(x-*)<"-i)-4°'^>i ^ \x*r' + \x2r-' ^ 21x2^-^ ^ 2 

\xi-X2\ ~ |a;*-a;2| |a;*| + |x2| ~ \x2\ ^ 

which is the inequality (5.6). 

We now turn to the inequality (5.7). When xi,X2 > or Xi,X2 < 0, it follows from a 
sharper bound — 03 1 < |ai — 02]^, valid for ai, a2 > andp G (0, 1). When > 0, a;2 < 
or xi < 0, X2 > 0, we have 

~ (x2)"-i| = Ixir^i + 1x21"-^ < 2(1x11 + I.T2I)"-' = 2|a;i - X2\''-\ 

Finally, the inequality (5.8) follows from |a:|" + \y\'^ ~ 2\xy\"/^ = _ |y|"/2|2 < + 

y|2"/2 = + and \x + y\" = |x + y|2"/2 = (|2:|2 + |y|2 + 2|xy|)"/2 < |a;|" + |?/|" + 2|a;y|"/2 
since < q;/2 < 1. □ 
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6. Application to central limit theorems 



The following result now follows easily from the covariance bounds in Theorems 2.1- 
2.4 and the central limit theorems for bounded functionals of SaS moving averages 
considered in Pipiras and Taqqu [16] and Hsing [11]. 

Theorem 6.1. Let a G (1,2) and {^n}„>o be an SaS moving average defined by (2.3) 
with a function a £ L"(R, dx) , a{x) = for x <Q. Suppose that K is a function satisfying 
Assumptions (K1)~(K3). // 



oo / /.m \ 1/2 oo 

J2( Hx^dx] <oo, X!Ko'^"]i <°°' 



(6.1) 



then 

N 



N-'/'Sn = ^-1/2 ^(i^(e„) - (C„)) ^ AA(0, a^) (6.2) 



n=l 

as N —> OO, where 

oo 

a' - lim Var(5jv) = Var(A'(eo)) + 2 V Cov(i^(eo), K{^n))- (6.3) 

N—^oo ^ — ^ 

n=l 

The limit in (6.3) exists and is finite and the series in (6.3) converges absolutely. 

Remark 6.1. The assumption of Theorem 6.1 that moving averages are causal (i.e., 
a(x) = for a; < 0) is also present in the available central limit theorems for bounded 
functions K (see Hsing [11], Pipiras and Taqqu [16]). 

The easiest way to verify condition (6.1) is to use, when possible, the stronger condition 
stated in the following lemma. Its proof is elementary and is given below. We use the 
notation [Co,Cn]7,5 = is,W{-x)\'^\a{n - x)\^ dx. 

Lemma 6.1. Let 7 > 0, 5>0 be such that 7 + (5 = a with a > 0. Also, let a G L"(R, dx). 
Then, 

00 / „m \ min{7,5}/a oo 

J2 i HxTdx) <oo =^ ^ [eo,Cn]7,5<oo. (6.4) 



— 00 ^ ^ n— — 00 



In particular, when a G (1,2), we have 

(a-l)/a 



°o y „jn >, (a-Lj/a oo 

( |a(x)|"dx <^ =^ ^[eo,Cn]i<oo, (6.5) 
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oo ^ .m \ 1/2 oo 

/ HxTdx) <c» ^ ^Ko,Cn]2<oo. (6.6) 

m=-oo\"'™-l ^ n=l 

Proof. By setting / = J2'^^i[io,£,n]'f.s, the result (6.4) follows from Holder's inequality, 

1= \a{z)\^\a{n + z)\^dz 

< ( l«Wrdz / Hn + zTdz 



E(/ i«wrd.) E(/ 



\a{z)\°'dz 

n+r/i— 1 



making the substitution n n — m in the last summation. □ 

Example 6. 1 . Let a G (1,2) and a : M ^ R be a bounded function such that a{x) = for 
a; < and |a(2;)| < Cx~p for x > and some p> 1/a. Then, for m>l and some constant 
C, 

\a{x)\°'dx<Cm-P"'. 

1 

Since < (a — l)/a < 1/2, the conditions on the left-hand sides of (6.6) and (6.5) hold 
as long as X^m^i < oo or, equivalently, p> l/(a — 1). Then, by (6.5), the 

function a satisfies condition (6.1) of Theorem 6.1 for p > l/(a — 1). 

Proof of Theorem 6.1. Suppose, for simplicity, that p=l in Assumption (K3) and 
set K(,{x) = A'(a;)l{|j.|>(,} and Ki^i,{x) = A'(x)l||3._^j|<i/f,|. Then, the function Lb{x) = 
K{x) — Ki,(x) ~ Ki_i,(x) is bounded. Denote the partial sum in (6.2) by Sb.N when K is 
replaced by Lb- To show (6.2), by Theorem 4.2 of Billingsley [6], it is enough to prove 
that: 

(i) N-^/'^Sb.N ^ A/'(0, as AT ^ oo; 

(ii) cr^ cr^ as 6 ^ oo; 

(iii) limsup5^g^ limsup^^^^ N^^ Yai{Sisi — Sb.N) = 0. 

Part (i) follows from Theorem 2.1 of Pipiras and Taqqu [16] since the first condition 
in (6.1) holds and the function Lb is bounded. Moreover, 

cr^= lim N-^YaT{Sb,N). (6.7) 



Part (ii) will follow from (iii). Indeed, since 



\VaY'^^iSb„N) - y&r'/^iSb,,N)\ <y&i-'^^{Sb„N - SN)+Vs^r'/\Sb,,N - Sn), 
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if (iii) holds, then ab converges to some a by the Cauchy criterion. To prove (iii), observe 
that 

N N 
n—l n—1 

and hence 

+ 2iV-i Var|^^^(A'i,b(e„) - Sifi,b(e„))^ 



= 2 Var(/vb(eo)) +4^1-]^ Cov{Kb{^o), Kbi^n)) 

+ 2Var(iCi,b(^o)) + 4 ^ - ^) Cov(ifi,b($o), ifi,fa(Cn)), 

n—l ^ 

which imphes that 

oo 

hmsupTV-i Var(5A, - Sb.N) < 2 Var(/fb(fo)) + 4 V | Gov{Kb{(o), Kb{U)\ 

OO 

+ 2 Var(A'i,b(eo)) + 4^ | Cov(Xi,f,(eo), Xi,6(e„))|. 

ri=l 

One can bound these sums by using Theorems 2.2 and 2.1 and Proposition 2.2. Sums 
involving the bounds [Co,?n]i converge by the second relation in (6.1) and the sum in- 
volving the bound [CojCn]2 converges by the first relation in (6.1) and Lemma 6.1 above. 
Part (iii) then follows by letting b^ oo. 

To show that cr^ is the asymptotic variance of N^^^^Sm in (6.3), first observe that 

Vari/2(^^) _ < |iV-i/2 Vari/2(^^ _ Sb^N)\ + \N-'^^ Ya^^^iSb^N) ~ <j\. 

Then, the first relation in (6.3) follows by taking limsup(,_^o2 limsup^^oQ in the above 
inequality and using (6.7) and parts (ii) and (iii) above. Finally, we need to show that 
the series in (6.3) converges absolutely. By Theorem 2.4 and Proposition 2.2, we get, for 
sufficiently large n, 

|Cov(/C(^o), A'(C„))| < C([eo,?„]i + Ko,^„]2). 

The absolute convergence then follows by using, as before, the relations (6.1) and 
Lemma 6.1. □ 
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Corollary 6.1. Theorem 6.1 is also valid for SaS moving averages {S,n}n>a defined in 
(2.3) by using function a such that a[x) = when x < Xo, for some xo < 0. 

Proof. Observe that the kernel a{n — x) in (2.3) can be replaced by a{n ~ {x — x^)) =: 
a{n — x) without changing the distribution of {^n}n>o- The function a is such that 
a{x) = for a; < and, if the function a satisfies (6.1), then so does the function a. Then, 
by applying Theorem 6.1 to the moving average defined via the function a, we conclude 
that (6.2) and (6.3) hold. □ 

Remark 6.2. A related paper of Wu [23] also contains a central limit theorem for 
infinite variance causal moving averages (not necessarily stable) of the "discrete" form 
Xn = ^'kLo^k£n-k and possibly unbounded functions K; see Theorem 1(a) in Wu [23]. 
The method of proof in Wu [23] is based on general central limit theorems for Markov 
chains developed by Michael Woodroofe and is quite different from the one presented 
here. The following observations shed light on the relationship between the result of Wu 
[23] and Theorem 6.1 above. 

First, there are unbounded functions K that satisfy the assumptions of Theorem 6.1 
but not those of Theorem 1(a) of Wu [23]. An example is the function K{x) ~ logj |a;| — 
E'logj 1^0 1, of particular interest in the wavelet-based application discussed below. The- 
orem 1(a) in Wu [23] involves the function Lk„{x), where Kn{w) = EK(w + Xn,i), 
Xn,i = aa£n + ai£„-i H h a„_iei and 

Lf[x) ~ sup 



y:\y~x\<l \y ^ 

is a local Lipschitz constant. For the logarithm function K above, 

Lk {x) > - — E = = oo 

ln2 \x + Xn,i\ 

and hence Theorem 1(a) of Wu [23] cannot be applied. 

Second, consider, for example, the power function K{x) = \x\^ — E\Xo\^ with /? < 
a/2, also of interest in wavelet-based application. Theorem 1(a) of Wu [23] applies to 
this power function under suitable conditions. We next examine how these conditions 
compare to those in Theorem 6.1 above. First, we can verify that the condition on Lk^ 
in Wu [23] is satisfied when (/3 — 1)2(7/ (? ~ 1) > ~1 or l/g < 2/3 — 1. Observe that since 
Wu requires g > 1, we obtain 2/3 — 1 > and hence a > 1, as in Theorem 6.1 above. 
Second, the condition on Kn in Wu [23] can be seen to be verified when (32q/ {q — l)<a 
or 1/q < {a — 2/3) /a. Comparing the conditions on Lk„ and K„, observe that they 
are verified if 1/q < 2/3 — 1 when 1/2 < jS < a/{a + 1), and if 1/q < {a — 2/3)/a when 
a/ia+1) <P <a/2. 

Now, suppose that the moving average coefficients a„ satisfy |a„| Cn^^ as ri — > cxd, 
where p > 1/a. The condition on a„ in Wu [23] requires that J2n Wnl"^^'^ < oo or 1< 
pa/2q. From the bounds on 1/q above, it is therefore necessary to have 1 < pa{2(3 — l)/2 
when 1/2 < /3 < a/{a + 1), and 1< p{a - 2/3)/2 when a/{a + 1) < P < a/2. On the 
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other hand, by Example 6.1 above, the conditions on a„ of Theorem 6.1 are satisfied 
when 1 < p{a - 1). When a/(a + 1) < /3 < a/2, observe that a — 2(3 < 2{a - 1). Hence, 
Theorem 6.1 is stronger than Theorem 1(a) in Wu [23]. When 1/2 < (3 < a/{a + 1), 
observe that a(2f3 — 1) < 2{a — 1) and hence Theorem 6.1 remains stronger. Combining 
these observations. Theorem 6.1 is stronger than Theorem 1(a) of Wu [23] for the above 
power function. Presently, we are not aware of any unbounded functions where Theorem 
1(a) of Wu [23] performs better. 

Theorem 6.1 has the following multivariate extension which is used in the next section. 
The proof of the extension is analogous to that of Theorem 6.1 and is omitted, except for a 
supplementary result on the form of a limit covariance. Let cj > be positive real numbers 
for j = 1, . . . , J. Consider the SaS moving average sequences {£,j,n}n>o, j = ^, ■ ■ ■ , J , 
given by 



S,j,n= / aj{n — Cjx)M{dx), (6.8) 
Jr 

where Oj € L"(M,da;) and aj{x) = for a; < and M is an SaS random measure with 

the Lebesgue control measure on R. Fix nj, j ~ 1, . . . , J, and let Nj be positive integers 
such that, as N —>■ oo, 

N 

Nj . (6.9) 



Theorem 6.2. Let a G (1, 2) and {£,j,n}n>o, j ~ ■ ■ i J ^ be SaS moving averages de- 
fined by (6.8) with aj such that aj{x) = when x < xq, for some fixed xo- Suppose that, 
for each j = I, . . . , J , the function Kj and the kernel aj in (6.8) satisfy conditions (Kl)- 
(K3) and (6.1), respectively. Then, asN^oo, 

(^7'^'%~.)/=i = NT^^^Y.^K,{i,^^)^EK,{i,,;)) Smo,a), (6.10) 

\ n=l / j^i 

where cr — {a'jk)j,k=i,...,j with 

(Tjk = hm E—:^^—, (6.11) 

N—too _/V N 
j k 

which exists and is finite. 

If, in addition, cj ~ 2~^ , Uj ~ 2^ and, for all k > j, X]^i[?i,"i?fc,o]i < oo, then the 
asymptotic covariance (6.11) is given by 

oo 

a^, = 2^i-k)/2 CoyiK,i^,,,,),Kk{^kfi)), k>j, (6.12) 

n— — OO 



where the series in (6.12) converges absolutely. 
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Proof. We prove the series representation (6.12). Recall that we now suppose Cj = 
and rij = 2K We shall also assume, for simplicity, that n^N^ = N in (6.9). (The general 
case can be proven in a similar way.) Then, Nj = 2^^^Nk and wc have, for k > j, 

m — 1 712 — 1 

= E E Cov(A-j(e,-„^_2.-,„J,iffc(a.o)) 

m — 1 n2 — 1 

= J2 J2 Cov(ii:j(^j-2'=-Jni-p-2'=-Jn2)'-f^fcfe.o)) 
p—0 ni — 1 712 — 1 

2'=-3-l Nk 

= E E (^fe-|"l)Cov(if,(e,,2'=-.„-p),J^fc(6.o)). 

p=0 n=-Nk 

After dividing by N^^nI^^ = 2^''~^'>/'^ Nk and letting iV -> c», wc then obtain by (6.11), 
that 

a,fc = 2(-'"-'=)/2 J2 E Cov(A',(C,-2^^.„_p),i^fc(&,o)) 

p—0 n— — oo 

oo 

= 2(^-k)/2 J2 Cov(A-,(0,„),/u(a,o)), 

n— — oo 

provided the scries is absolutely convergent. This can be established by first using a result 
for (,j,n, analogous to Proposition 2.2 and Theorem 2.4, to conclude that, for sufficiently 
large n, 

|Cov(K,(C,,„),J^fe(&,0))l<C(K,,n,6.0]l + [Cj,n,?M]2). (6.13) 

When k = j, the convergence of the series with the terms on the right-hand side of (6.13) 
follows from the fact that the aj's satisfy (6.1). When k > j, we use the assumption 
EnKi:"'^'=-o]i < oo. together with the observation that the first relation in (6.1) with 
a = aj and a ^ implies ^„[Cj,n7^fc.o]2 < oo, which can be proven as was (6.6) in 
Lemma 6.1. □ 



7. Asymptotic normality of wavelet-based estimators 

We here apply Theorem 6.2 to prove the asymptotic normality of wavelet-based estima- 
tors of the self-similarity parameter in linear fractional stable motion (LFSM) . LFSM is 
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an SaS self-similar stationary-increments process with the integral representation 

^it)^ [ {it-u)l-i-u)l}M{du), teR, (7.1) 

where 

K, = H -1/a, (7.2) 

H € (0, 1) is the self-similarity parameter, a S (0,2) is the stability parameter and M is 
an SaS random measure with the Lebesgue control measure on K. Self-similarity means 
that for all c > 0, the processes X{ct) and X{t) have the same finite-dimensional 
distributions. LFSM is an infinite variance counterpart of fractional Brownian motion, 
which is the only Gaussian self-similar process with stationary increments. It is often 
taken as a representative process for self-similar processes with stationary increments 
having infinite variance. (The more general double-sided LFSM cannot be considered for 
Theorem 7.1 below because of the causality assumption discussed in Remark 6.1.) For 
more information on LFSM, see Section 7 of Samorodnitsky and Taqqu [19]. 

Wavelets have already proven useful for estimating the self-similarity parameter of 
fractional Brownian motion and a related long-memory parameter in finite variance long- 
memory time scries (Vcitch and Abry [22], Bardct [5] and Abry, Flandrin, Taqqu and 
Veitch [2]). This work was motivated to a great extent by applications to data traffic 
in communication networks. Following this line of research, several authors have sug- 
gested using wavelets to estimate the self-similarity parameter H in LFSM as well. From 
a mathematical point of view, it is interesting to see to what extent existing wavelet 
methodology applies to this extended class of processes. There is indeed evidence that 
some teletrafhc data deviates from Gaussianity and has heavy tail, characteristic of pro- 
cesses with infinite variance. 

The discrete wavelet transform of LFSM (or of other deterministic or stochastic func- 
tions) is a sequence of discrete wavelet transform coefficients {dj^k}j,kez, defined by 

dj,k^ ( X{t)i^,,k{t)dt^ [ X{t)2-''^i:{2-H-k)At. (7.3) 

JR 

Here, : R i— > R is the so-called wavelet function which has Q > 1 zero moments, that is, 

I tp{t)t"'dt^O, m^O,...,Q-l and / V(0^'^ dt ^ 0. (7.4) 
Jr Jr 

It was shown by Delbeke and Abry [8, 9] and Pesquet-Popescu [15] that the discrete 
wavelet transform coefficients of LFSM are well defined for a bounded, compactly sup- 
ported wavelet ip when H — \/a> —1 and that the following Fubini-type result holds for 
j, /c e Z: 

d,,k^ I ( I {t~u)l4,j^k{t)dt\M{Au)^ f 2^(-+'/^^h{k~2-^u)M{du), (7.5) 

JR \JR J JR 
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where 

h{u)^ J {s + u)ltl;{s)ds, u€R. (7.6) 

Moreover, it is easy to see that for each j G Z, the sequence {dj^k}kei, is stationary, that 
is, for all I G Z, {dj j^^i} =rf {dj k}, and that the following scaling relation holds: 

K,},,z = {2^^^+^/''do.fcK.gz. (7.7) 

Let us also note that, in practice, the wavelet coefRcients are computed (rather, approx- 
imated) by using fast pyramidal Mallat-type algorithms (see, e.g., Mallat [14]). 
A wavelet-based estimator of H can be defined by 

1 1 
H^J2'^jj^J2^og,\d,,„\~-. (7.8) 

j n — 1 

The summation here is over some finite number of consecutive indices j, called octaves 
or scales, and each integer Nj corresponds to the number of available wavelet coefficients 
of LFSM over the time interval [0,-/V] at scale j. Roughly speaking, we have Nj = 2~^N 
(up to border effects) and thus Nj — Nj{N) oo as oo. Finally, the Wj's satisfy 

the relations 

Y.w,=Q, ^j-i«, = l (7.9) 

j j 

and can be viewed as weights for a linear least square estimation of a slope in {i,Yj), 
where Yj is the term multiplying Wj in (7.8). While the wavelet coefficients dj^k have 
infinite variance, log2|rfj^A;| has finite variance, as do the estimators H. Theorem 6.2 
implies that this estimator is asymptotically normal. 

Theorem 7.1. Let H be the self-similarity parameter estimator in LFSM defined by 
(7.8) for bounded and compactly supported wavelets. Suppose that a G (1,2) and that 

Q~H> ^ (7.10) 
cx.{a — 1) 



Then, 



as N ^ oo, where 



tith 



N^'^{H -H)^U{Q,(J^) (7.11) 



= 2(^-'=)/2 Cov(log2|d,-„|,log2|dfc,o|) (7.12) 
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if k > j , and CFj^k = <^k.j if j > k. We also have 

^ lim NE{H-Hf. (7.13) 

Remark 7.1. Because convergence in (7.11) is to a normal law, a rate of iV^/^ is ex- 
pected. In Kokoszka and Taqqu [13], the Whittle method is used to estimate the long- 
range dependence parameter in stable FARIMA time series. A faster rate of essentially 
N^/°' is obtained, but also a limit which has infinite variance and is hence more spread 
out. 

Proof of Theorem 7.1. Without loss of generality, we suppose that the sum in (7.8) 
is over octaves j ' = 1, . . . , J for some fixed J. By using (7.7), we can write 

1 

7Vi/2(^_^) =^^^.2.72__^(log^ |d,,„| -i?log2|d,,„|) =wY*, (7.14) 

where t stands for "transposed," = [i^/'^wi,. . . , 2'^/'^wj) and Y = (F(l, iVi), . . . , Y{J, Nj)), 
with 

Y{j, N, )^NJ^'^Y. (log2 \dj,n\~E log2 I d,- „ I ) . 

Therefore, to prove that H is asymptotically normal, it is enough to establish the asymp- 
totic normality of the vector Y. Observe that, by (7.5), this vector can be expressed as 
that in (6.10) and (6.8) with ^ and aj{x) = 23'-'^+^^^^h{x), where j = 1,...,J. 
Moreover, since the wavelet ip has compact support, there exists an xq G M such that 
h{x) = for X < Xq and thus aj{x) = for x < xq as well. Then, by Theorem 6.2, we get 
asymptotic normality, provided the function log2 |a;| satisfies (K1)-(K3) and the function 
a{x) = h{x) satisfies (6.1). It is easy to verify that K{x) ~ log2 |a;| satisfies (K1)-(K3). 
The function h satisfies (6.1) by using Lemma 7.1 below and Example 6.1, provided that 
Q — K> l/{a — 1) or, equivalently, Q ~ H > l/{a{a — 1)). 

Let us now show that can be expressed as in the theorem. By using (7.14) and 
Theorem 6.2, it is enough to show that under the assumption (7.10), '^'^^i[dj^m dkfi]i < 
oo for all k > j. This can be established as in Example 6.1 by using aj{x) = 2^^^^-^/^^h{x), 
Lemma 7.1 below and the implication (6.5). □ 

The following result was used in the proof of Theorem 7.1 above. 

Lemma 7.1. Let h be the function defined by (7.6), where the wavelet ip has compact 
support, is bounded and has Q>1 zero moments. Then, for sufficiently large u, \h{u) \ < 
Cu'^~'^ , where C is some constant. 
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Proof. By applying Taylor's formula to the function fu{s) = (s + u)'^, which is infinitely 
dilFerentiable on any interval s S [— M, M], M > 0, for sufficiently large u, we have that 



- 1)---(k-Q + 1) q q 



where sq £ [— M, A/]. If the wavelet has compact support, we then deduce the result 
by using (7.4). □ 

Note that the greater the value of Q in Lemma 7.1, the faster \h{u)\ decreases as u 
grows. Thus, in view of (7.5), when Q increases, the sequence {dj^k}kei. becomes almost 
independent. This is one of the main advantages of working with the wavelet coefficients 
rather than with LFSM itself A more comprehensive and applied study of wavelet-based 
estimators of the self-similarity parameter in LFSM (and further results) can be found 
in Abry, Delbeke and Flandrin [1] , Abry, Pesquet-Popescu and Taqqu [4] , Abry, Fladrin, 
Taqqu and Veitch [3], Stoev, Pipiras and Taqqu [21] and Stoev and Taqqu [20]. 

Remark 7.2. An alternative estimator for the self-similarity parameter in LFSM can 
be defined as 



(1 \ 1 

■I n— 1 / 



F* = ->^u;,log2|— > Jd,,„|'^)--, -l<p<a/2. (7.15) 
The asymptotic normality result for H* can be established as in Theorem 7.1 with 



n— — oo 



if k > j , and dj.fc = a^j if j > k. However, because of the presence of the function log2 in 
(7.15), we cannot deduce that the relation (7.13) holds, namely that cr^ is the asymptotic 
(normalized) variance of H*. 

Remark 7.3. Another interesting question not addressed here is that of joint estimation 
of H and a. We plan to pursue this question in future work. 
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